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Abstract 

We have studied extracting iV^bAcbl by calculating the ratios B(B° — > 
D^~(Tr + ,p + ))/B(B° -> D ( S * ] ~D^ + ) including penguin effects within the 
factorization assumption. The ratios involving B° — > D~D + mode have con- 
siderable penguin corrections (~ 15% at the amplitude level), but those in- 
volving ~E° — > D~D* + mode have relatively small penguin corrections. On 
the other hand, the 5° — > D~D + mode has smaller form-factor dependance. 
Therefore, these ratios complement each other in measuring Vub/Vrfi- Tbe 
theoretical uncertainty from the hadronic form factors in our method is at 
the level of 15%, which is comparable to the model-dependence uncertainty of 
about 20% in the measurement of (V^b/V^bl from the exclusive semileptonic 
B decays. Using the newest upper limit on B — > D s tt decay from CLEO, 
our method sets an upper limit (V^/Kbl < 0.13 which is very close to the 
measured values from the semileptonic B decays. We also discuss the possible 
breaking of factorization assumption. 



I. INTRODUCTION 



A precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [||] is 
one of the key issues in the study of B mesons and S-factory experiments 0. For instance, 
the value of V u \> being non-zero is a necessary condition for CP-violation to occur in the 
SM; otherwise, we have to seek for new physics explanation for the observed CP violation 
in the Kl — > tttt decays ||. For a stringent test of the SM and a search for new physics, 
it is important to make a precise measurement of the modulus |Kib| from the decays of B 
mesons. 

Theoretical and experimental studies on have been mostly focused on the semi- 
leptonic B meson decays. Observations of inclusive and exclusive semileptonic b — > u tran- 
sitions by the CLEO f|,|] and ARGUS H experiments confirm that V n h is indeed nonzero, 
but these measurements suffered from large uncertainty due to model-dependence in extract- 
ing the value of V n b- For instance, determination of V^b from the inclusive b — > uiv using 
the leptonic end-point momentum involves large uncertainty in determining the fraction of 
partial decay width for Eg being larger than some cut value. In particular, the dependences 
of the lepton energy spectrum on perturbative and non-perturbative QCD corrections as 
well as on the unavoidable model-specific parameters are the strongest at the end-point re- 
gion, which makes the determination of |Kib/Kb| very difficult by this method. There have 
been other suggestions for avoiding these difficulties in studying the inclusive semileptonic 
b — > u decays, for example, by using the invariant mass of the hadronic system |7]||] recoiling 
against £z/, or the invariant lepton mass |j| . Experimental studies have been made using the 
invariant mass of the hadronic system to separate b — > u from b — > c decays, and measure 



|Kib/Kb| IH3- By considering the exclusive semileptonic decay modes such as B — > iriv 
and p£i>, we may reduce the dependence on the Eg spectrum, but we are confronted with 
different aspects of theoretical uncertainties. In this case, determination of V^b becomes very 
sensitive to the fraction of such exclusive decays with respect to the inclusive b ulv and 
it has large uncertainty coming from hadronic form factors. 
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Although traditional difficulties with the understanding of non-leptonic weak decays have 
prevented their use in determination of CKM elements, the possibility of measuring |V^b| 
via non-leptonic decays of B mesons to exclusive two meson final states |TT| , |T2"11 has been 
theoretically explored. To avoid the theoretical difficulties of non-spectator decay diagrams, 
only those final states have to be chosen in which no quark and anti-quark (qq) pair has the 
same flavor W2 . 

Within the factorization approximation^ and after considering the final state interactions, 
exclusive two body decay modes of B mesons would certainly be worthy of full investigation. 



In Ref. j llj, it was pointed out that for the extraction of |Kib/Kb| from nonleptonic B 
meson decay data, study of the ratio B(B° — > D+n~)/B(B° — > D+D~) is useful. Since 
the final states 7r~.D+ (D~D+) consist of a single isospin component / = 1 (/ = 1/2), 
the decay amplitudes are independent of the phase shift caused by elastic rescatterings in 
final states. Furthermore, the decay mode B° — > D+tt~ is caused by only one diagram, 
the b — > u tree transition with external W emission and so completely independent of the 
penguin-type interaction. On the other hand, in the decay mode B° — > D+D~, although 
the dominant contribution is from the b —>■ c tree diagram, b —>■ s penguin diagrams also 
contribute. In Ref. [nj the authors neglected the penguin effect by expecting its size to be 
small compared to that of the tree diagram. However, we find that the theoretical estimate 
of the penguin correction is not small enough to be simply neglected for B° — > D+D~ decay 
mode. On the other hand, we find that B° — > D+D*~ decay mode receives much smaller 
penguin corrections compared to B° — > DfD~ decay mode. In this work we investigate 
more thoroughly the possibility of extracting iKib/Kbl from study of the nonleptonic decay 



ratio B(B° -> £><*>- (tt+ p + ))/B(B° -> £>W-£>M + ). 



For possible breaking of the factorization approach, we give a short discussion in Section II. 



2 



II. THEORY ON B° -► (vr+,p+) AND 5° -► D { g ] ~ £)(*)+ 

AND NUMERICAL ANALYSIS 



While only tree diagrams contribute to 5 — > D~n + decay in the SM, the decays S — > 
_Dj £)(*)+ receive penguin contributions as well. The relevant A.B = 1 effective Hamiltonian 
has the form; 



v/2 V 



cs 



Cl (/x)0?(/i) + c 2 (/i)0 2 V) - ^e^mm + h - c -> (!) 



i=3 



where 0± i2 represent QCD corrected tree-level operators and 3 _ 6 (O 7 _i ) the QCD (elec- 
troweak) penguin operators, which are defined as 

Of = q-jf.LbsYLc, 0\ = qa^Lbps^Lca, 

03(5) = s^Lbcry^^R)^ 4 ( 6 ) = s fl li U)pcp'fL{R)c a , (2) 
07(9) = s^Lbcry^R^Ljc, O 8(10 ) = s a -i^Lbpc^R(L)c a , 

with i?(L) = 1 ± 7s and g = u (c) for 5 -> D s tt (5 -> D S D^>) decays. 

In the factorization approximation, the decay amplitudes of our interest are expressed as 

A (B° _> £>(*>-(*+,+)) = ^K^ c ;a 1 ( J D;|s7^c|0)(7r + ,p + | M7 ,L6|5°), (3) 



4(5° - = % (k 6 K>i - W «4 + a 10 + 2(a 6 + a 8 



x( j D;|s7^c|0)( j D + |c 7m L6| j B ), (4) 



A(B° -> D~D* + ) = % ( - Kft^ a 4 + a 10 - 2(a 6 + a 8 ) 



m 



x( J D;|s 7 ^c|0)( J D* + |c7^|fi°), (5) 

A (B° _> = ^ {W c >i - ^(04 + a 10 )} (£>r |s 7 ^c|0)( J D + |c7 M L6|B°) (6) 

Here a/s represent effective parameters defined as 



where cf 1 are the renormalization scheme and scale independent effective Wilson coefficients, 
and N* s is the so-called effective color number, which is supposed to include non-factorization 
effects as well as color suppression effect, and can thus be considered a free parameter. 
Using VtbV t * = — V c bV* s , one can cast the amplitudes in more compact forms 

A(B° -> #(*)-£>(*)+) = ^fv cb V c *M B -> D^W^)(D^-\sYLc\0}(D^ + \c lfl Lb\B°} } (8) 

V2 



where 



a 1 (B^DD s )=a 1 [l + ^^ + 2 



0-4 + OlO , o a 6 + «8 



5i(S -> D*D S ) = ai 1 + u - 2 



ai ~ ai (m& - m c )(m c + m s ) 
a 4 + ai a 6 + a 8 mf> 



ai ai (mj + m c )(m c + m s ) 



a^B- £>£>;) = ai (l + ^+^V (9) 



Using the numerical values of Oj's in Ref. [DJ, the effective parameters o4 defined above are 
related to a± by 

|a x (B -> £)£) s )| = 0.847a!, 
|5i(S -> £>*£> s )| = 1.037a!, 

|oi(S->iPD;)| = 0.962a!. (10) 

From the above relations one can see that, at the amplitude level, the penguin contributions 
to 5° — ► D~D* + decay (3.7%) are much smaller than those for B° — > D~D + mode (15.3%). 
Actually the penguin effects on B° — > Dj D + decay are not small enough to be simply 
neglected. As mentioned in the Introduction, the penguin effects are neg lectedinRef. Q. 
We note that for B° — ► D~D* + decay mode the penguin contribution can be neglected. This 
difference of penguin contributions to the similar modes 5° — ► DjD + and — > D~D* + is 
due to the different chiral structure of the final states. B — > D* transitions occur through 
axial vector currents, while B — > D through vector currents. 
Then, the ratios 



B(F-+D 7 (n+,p+) ) 
B(B^ -> D-D&+) 



^•(7T,p)/D(*) = — " ~ (11' 



and 7^ /I3 = -^zjj 7 12 

£(F -> £>;-£>+) 
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are given as 



ft 



n/D 



Kb 2 



V, 



cb 



a x {B^DD s ) (m 2 B -m 2 D ) 2 L C D , 



T?B-n 
^0 



(13) 



ft 



p/D 



ub 



cb 



m 



P, 



p3 ' 



2 \2 „D 



(14) 



ft 



7r/D* 



ub 



v; 



cb 



ai 



— m; 



m. 



2\2 



L>*3 



m 



2A 



BD* 



m 



(15) 



ft 



V, 



ub 



cb 



^(5 -> £>*£> s 



D* 



( < P ^ls 



m 



Ds 



(16) 



ft-vr/ 



ub 



Vr 



cb 



ai 



di(B -> DD*)J \tf \F? D {m%). 



(17) 



where p;f is the cm. momentum of the decay particle X. Here the form factors follow the 



following parameterization WA 



[P\p')\VAP{p))=[p,+p^ 

2 



m 2 p — m 2 P , 



m P — m,p/ 



(Vy,e)|V M |P(p)) 
(V(p',e)\A,\P(p)) 



VTip + VTly 



e^e*»p a p ,(3 V(q 2 



[m P + m v )e /i Ai(g ) - 



mp + my 



2m v e -f qfM [A 3 (q 2 ) - A Q (q 2 )} 



where q = p-p', F 1 (0) = F (0), A 3 (0) = A o (0), 

2 nip + my 



Ml* 



2m, 



m P — my 



2m 



A 2 (q 2 



v 



and P, V denote the pseudoscalar and vector mesons, respectively. 

As can be inferred from the explicit expressions in Eqs. (]T3|) - Ql7D , these ratios do not 
have much dependence on values of a^. Especially, for 5° — > DjD* + decay, the penguin 
contributions add up destructively and so the dependence on a, is almost cancelled out in 
the ratio. However, explicit calculations for the ratios of branching fractions depend strongly 



on the form factors. In the following analysis, we consider five models for the form factors of 



B — > 7r transitions: two quark-model approaches (the Bauer- Stech-Wirbel (BSW) model |14 
and Melikhov/Beyer pf), light-cone sum rules (LCSR pj), lattice QCD (UKQCD | 



and relativistic light-front (LF) quark model [18|]. And for B — ► transitions, we adopt 
BSW, Melikhov/Stech |TjJ, and LF models. In Table 1, we list explicit numerical values of 



m 



K 



tt/D 



B(B 



D 



form factors evaluated at q 2 

B(B° -> D-7i 



(*) 



. Then we get the theoretical predictions 



D-D- 



0.424 



Kb 2 



V, 



cb 



m 



1 2 



0.319 



0.740 



1 2 



[0.424 ±0.041] 



ub 



K 



cb 



(19) 



ft 



p/D 



B(W -» D.y 



0.443 



v; 



ub 



K 



cb 



-, 2 



n? 



0.398 



[0.443 ± 0.063] 



V, 



ub 



Kb 



0.740 



-^0 l m D s 



(20) 



B(B° 



B(B° -> D~D* 



0.459 



ub 



cb 



/•',;•" ( ;// 



[0.459 ± 0.076] 



0.319 



Kb 2 



K 



cb 



0.793 



(21) 



-7? 

/<-p/D* = — =q- 



B(B -> 



0.480 



K 



ub 



m 



Kb 

[0.480 ± 0.094] 



0.398 

ub 



K 2 



Kb 



0.793 



AE D *{ml s ) 



(22) 



ft 



B{B -> J D*- £)+ ) 



0.456 



k 



ub 



K 



cb 



0.367 



[0.456 ± 0.038] 



ub 



v; 



cb 



2 r 



0.817 



l 2 



(23) 



where the errors are originated only from the dependence on the hadronic form-factors. 
Considering the current experimental results [pOlpTl 



B(B -> DJtt + ) < 5.1 x 10" 



B(B^ 
B(B° 

23(5° 



£>7p + ) < 7.0 x 10~ 4 , 
D-D + ) = (8.0 ±3.0) x 10" 3 , 
D-D* + ) = (9.6 ±3.4) x 10~ 3 , 
D*-D + ) = (1.0 ±0.5) x 10~ 2 , 
D*~7T + ) < 7.5 x 10~ 5 , 



(24) 



we then estimate 



V, 



ub 



v r 



cb 



< < 



0.13 


from 




0.12 


from 


TZtt/D* 


0.14 


from 




0.48 


from 




0.44 


from 


IZp/D* 



(25) 



We note that with the new preliminary upper limits on B — ► D s ii + from CLEO f21 |, the 
upper limit on |Kib/Kb| is already very close to the current estimate |Vub/Kb| = 0.09±0.025 
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Next we consider direct extraction of |V^b| from the B° — > D~n + decay rate. As men- 
tioned earlier, this mode is one of the cleanest processes without any penguin corrections. 
Main uncertainties come from the form factor Fq n and Wilson coefficients or the effective 
parameter a 1 . The decay rate is given as 



(1.065 x 10 



ri2 

2 li; " r ' 



l T/*|2 Pc 
ub v cs I o 2 



ml - K) 2 [a 1 f D .F**{ml 



B 



- 12 Wu h \ 2 



a i 



1.059 



D s 



0.240 GeV 



Bit 1^,2 
Ds 



0.319 



GeV. 



(26) 



In order to estimate theoretical uncertainty from the effective coefficient a±, we choose dif- 
ferent values for Nfr , the effective number of colors appearing in Eq. (]?]), Nfr = 2, 3 and 
oo, and we use for the renormalization scheme independent Wilson coefficients = 1.149 
and c| ff = —0.325 [[HJ. Considering the new preliminary upper limits on B° — ► Djn + in 
Eq. ([24]), we then estimate 
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I Kb I < (0.450 ±0.020 ±0.021) x 10" 



(27) 



where the first error is due to the effective coefficient a± and the second one is from the 
dependence on the form factor F^ n . 



Finally we note on the possible breaking |23[ of the factorization assumption, on which 
our previous results are based; we now generalize Eqs. (|i"3|) - ([D]) by including the factor 
for breaking of factorization, 

2 r F B-,ii 



£(£ L 



ub 



cb 



K„ 



i/j 



a i 

a, i 



i/j 



[m 



D 



(*), 



F B^J( m 2 



D 



(*)/ 



1 ± (F — B) 



i/j 



(28) 



where i = n,p and J = D,D*. Ki/j is the kinematic phase space factor, and a\ja\ is 
the ratio of effective Wilson coefficients, for which we used values estimated within factor- 
ization assumption. Because the decay — > D^~~i + is only through tree diagrams and 
IT -> DW-J+ is polluted by penguin diagrams, the ratios IZi/j would possibly be the 
best observables to measure the breaking of factorization, (F-B)^j. Instead of measuring 
iK&Ac&l as proposed, if we use the value of iKft/K&l measured from the semileptonic decays 
of B mesons, then we can systematically estimate the (F-B)j/j through experimental values 



of n 



i/j- 



III. DISCUSSIONS ON EXPERIMENTAL FEASIBILITY AND SUMMARY 

We have investigated the possibility of extracting | Kb/Kb | and |Vub| from non-leptonic ex- 
clusive decays of B meson into two meson final states. In particular, we calculated the ratios 
B(B° — > Z)M~(7r + , p + ))/B(jf — ► D^~D^ + ) including penguin effects in the factorization 
assumption. By taking the ratios, some model-dependence on the coefficients Oj and hadronic 
form-factors is reduced. We found that the B° — > D~D + mode has considerable penguin 
corrections (~ 15% at the amplitude level) which cannot be simply ignored as it was done 
in Ref. ||11|| . We also found that the B° — > D~D* + mode has very small penguin correc- 
tions. On the other hand, the 5° — > D~D + mode has smaller form-factor dependence than 
— > D~D* + . Therefore, these modes can complement each other in measuring V n b/Vcb, m 
conjunction with B — > D^n and B — > D s p. 
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We have shown that the relevant hadronic form factor uncertainties in our methods are 
typically at the level of 15% (Eqs. (|19|) - (|23|)) . On the other hand, the mo del- dependence 
uncertainties in the measurement of |Vub| from exclusive semileptonic B decays are currently 
at the level of 20% 0. Therefore, if we can contain the experimental uncertainty of our 
method within 15%, the method described in this paper becomes competitive with the 
semileptonic analyses. This implies that we have to determine the decay branching ratios of 
B — > D^*'tt, etc. with experimental uncertainties being less than 30%. 

Current experimental uncertainties for the branching ratios of B — > D S D* decays are at 



the level of approximately 30%. The recent CLEO measurements |24| of B(B — » D s D + ) = 
(8.4 ± 3.0) x 10~ 3 and B(B° -> D;D* + ) = (9.0 ± 2.7) x 10~ 3 are based on the event 
sample of (2.19 ± 0.04) x 10 6 BB pairs. With high-statistics event sample expected from 
the -B-factories, we anticipate reducing the uncertainties of these decays to below 10% level 
in the very near future. As for the decay B — > D^n and D s p, none of them have been 
experimentally measured yet. Recently, CLEO has presented the following preliminary upper 
limits: B(B -> D s n + ) < 8.9 x 10~ 5 and B{B -> D*ir + ) < 7.5 x 10~ 5 fgj], based on an 



event sample of 9.7 x 10 6 BB pairs. We expect to have larger amount of data from each B- 
factory experiment very soon. This, combined with much improved hadron identification and 
vertexing capabilities of the .B-factory experiments, will improve the sensitivity of B — > D^n 
and D s p searches down to B ~ 10~ 5 level and discovery of the B — > D^Ti(p) reactions might 
be possible in the near future. 

Using the current estimate of | V^/Kbl = 0.09 ± 0.025 p2]| , along with the world-average 



value of B(B° -> D~D+) = (8.0 ± 3.0) x lO" 3 |3, we take B(B U -> Djir+) = 2.7 x 10" 5 as 



-ftO 



the expected branching ratio. Based on this, we estimate experimental conditions to achieve 
30% statistical uncertainty in B(B° — > D~n + ) from the -B-factories. Consider we measure 
B — > D s n decay through D s — > (pir and <fi — > K + K~ . In terms of signal detection efficiency 
e and integrated luminosity £ (in fb _1 ), the expected number of reconstructed B — > D s 7r 
events are 

iV = e £ x 10 6 x 2.7 x 10" 5 x 0.018 



= 0.5 e C, 



where 0.018 is the decay branching ratio of the particular D s decay mode that we consider. 
Assuming e = 17% |25| we need C ~ 120 fb" 1 to obtain 10 signal events. We can further 
improve the experimental sensitivity by several factors, if we include other decay channels 
of D s and if we also analyze related modes such as B — > D*tt, B — > D^p, B — > Dy'u}, etc. 
In this case, with i2 ~ 50 fb _1 of data from the 5-factories, it may be possible to obtain a 
competitive, independent measurement of from exclusive non-leptonic B decays. 

There have been other studies of using nonleptonic B decays for measuring \ V u b/V c b\. For 
example, there were studies of utilizing fully inclusive nonleptonic b — > ucs' decays ]26| or 
exploiting semi-inclusive nonleptonic B decays B — > DfX u [^7|. The method suggested 
in Ref. is clean in theoretical calculation. But it is based on the assumption that 

we can experimentally separate, without introducing much model-dependent uncertainty, 
b — > ucs' decays from the dominant b — > ccs' decays which is two-orders-of-magnitude larger 
in size. With existing experimental capabilities of I?-factory experiments, this assumption is 
not justified. Ref. |2^J considers more carefully the experimental backgrounds, but ignores 
the significance of the continuum background which can be statistically dominant in the 
interesting range of 2.0 < po s < 2.5 GeV/c 
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TABLES 

Table 1. Numerical values of form factors at q 2 = m 2 w in various form-factor models 







F^{m 2 Dt ) 


F BD (m 2 Da ) 


Ff D (m 2 Dl ) 


A$ P (m 2 Ds ) 




BSW 


0.377 


0.395 


0.753 


0.776 


0.326 


0.690 


Beyer /Melikhov 


0.299 


0.357 






0.381 




LCSR 


0.322 


0.382 






0.457 




UKQCD 


0.301 


0.359 






0.469 




Melikhov/Stech 






0.722 


0.806 




0.818 


LF 


0.295 


0.349 


0.746 


0.869 


0.357 


0.872 
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